Уезтік гооіояії, 40(2): 99—113, 2006 
о Л. И. Францевич, 2006 


УДК 502.5:621.039.7 


СТАНДАРТИЗАЦИЯ ДАННЫХ ДЛЯ СОСТАВЛЕНИЯ 
КАРТЫ РАДИОАКТИВНОГО ЗАГРЯЗНЕНИЯ ЖИВОТНЫХ 


Л. И. Францевич 


Институт зоологии им. И. И. Шмальгаузена НАН Украины, 
ул. Б. Хмельницкого, 15, Киев, 01601 Украина 


Принято 3 марта 2006 


К 20-летию Чернобыльской катастрофы 


Стандартизация данных для составления карты радиоактивного загрязнения животных. Франце- 
вич Л. И. - На карту окрестностей Чернобыльской АЭС (85 х 65 км) выносятся данные о 
загрязнении '7С$ рыб, птиц, грызунов и копытных и данные по бета-активности раковин 
моллюсков как показателе загрязнения ?05г. Задача обработки данных — свести массив 
измерений к стандартной дате и к видам, принятым за стандарт. Установлены сезонные тренды 
радиоактивного загрязнения для кабана (,5и5 5сго/а) и косули ( Сарғеоіиѕ саргеоіиѕ), многолетний 
тренд для моллюсков и рыб, стабильные отношения накопления у разных видов, родов, 
экологических групп у моллюсков, рыб, копытных. Данные измерений пересчитаны на 
стандартный вид (для моллюсков — большой прудовик, Гутиаеа 5авпай5; для рыб — карп, 
Сурғіпиѕ сагріо; для копытных — кабан, 5и5 5сго/а) и стандартную дату 01.10.1996. Стандартиза- 
ция данных вдвое уменьшает первоначальную дисперсию измерений в данном местообитании. 
Составлена номограмма для вычисления процента образцов, превышающих данное значение 
загрязнения, при известных значениях среднего и коэффициента вариации. Установлены 
коэффициенты накопления и перехода. Обнаружена положительная корреляция загрязнения 
проб из соседних точек местности, тем большая, чем меньше расстояние между пробами. 
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аа 5іапдагіїгабоп юг Марріпє ої Кайіоасііуе Сопіатіпайоп оѓ Апітаї8. Егапбеуей Г. І. — Тһе 
тар ої уісіпібіеѕ ої Ше Срегпобу! Миеаг Ромег Ріапі (85 х 65 кт) Ферісіє даа оп сомапитаНоп 
оЁ ће ћѕһ, Ыгаѕ, гойепіѕ, апа ипешаѓеѕ \мИВ З"С8, аз ме! аз даа оп е Беѓа-асііуіѓу ої тоПиѕс вез 
іпаісаіпе сопіатіпайоп ув °5$г. Тһе роа! ої Чайа ргосеѕѕіпо \аз їо гедисе һе аггауѕ ої 
теаѕшетепіѕ їо ѕоте 51апдага ааѓе апа ѕѓапаага ѕресіеѕ. Зеуега! геп4$ Гог гайіоасіуе сопіатіпаїйоп 
Пауе Бееп геуеаІеа: Ве зеазопа! їгепаѕ іп һе Боаг (5и5 5сго/а) апа ће гое 4еег ( Саргеоійз саргеоіиѕ), 
Фе 1етрога! регеппіа! ігепа іп тоПиѕсѕ апа #ёѕһ, 5габіе гаііоѕ ої асситиіайоп іп аіегепі ѕресіеѕ, 
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ѕќапаага ѕресіеѕ (ропа зпаЙ Гутнаеа ѕіарпаііѕ атопе тоПазК$, сагр Сургіпиз сагріо атопе, ћѕћ, Боаг, 
Уиз 5сго/а атопе ипешаѓеѕ) апа даже 01 Осіобег 1996. Тһоѕе ргоседигез йітіпіѕћ (ће іпійа! уагіапсе 
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Кеу мога: гойепіѕ, ипецаѓеѕ, тоПиѕкѕ, Ыігаѕ, сагіовгарпу, асситиіайоп Гасіог, їгапѕѓег стог, 
гайіоесоІову, Срегпобу! Саѓаѕігорће. 


Постановка задачи 


Институт зоологии им. И. И. Шмальгаузена НАН Украины выполнял картирование радио- 
активного загрязнения животных как часть проекта «Атлас Чернобыльской зоны отчуждения», первый 
том которого издан в 1996 г. (Атлас..., 1996). Есть два подхода к картированию собранных данных: 
1) разместить на карте все точки сборов с указанием первичных радиометрических данных для всех 
видов животных (карта фактов) или 2) создать карту как инструмент прогноза загрязнения животных. 

Трудности картирования радиоэкологической ситуации возрастают по мере перехода от карт 
загрязнения почвы и водных объектов к картам загрязнения растений и животных. Животный мир 
Зоны отчуждения состоит из 400 видов позвоночных и нескольких тысяч видов беспозвоночных с 
различным образом жизни и биохимическими особенностями. На карте, однако, может быть отобра- 
жено пространственное распределение усредненных значений немногих величин. 
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Изменчивость радиоактивного загрязнения — объективная особенность животных. Так, загря- 
знение рыб !37С$ в пойменном озере Глубокое изменяется в широких пределах от 2 до 130 кБк/кг 
(июнь, 1995), в пруде-охладителе от 4,4 кБк/кг (карп, ноябрь, 1995) до 318 кБк/кг (судак, июнь, 
1987). В сборах рыжей лесной полевки из Янова найдены особи с загрязнением от 12 175 кБк/кг 
(05.09.91) до 265 кБк/кг (19.09.93). Пределы загрязнения для кабанов, добытых в Кулажине, 
составили 4334 кБк/кг (24.04.92) — 5,5 кБк/кг (03.11.93). 

Состав кормов меняется с возрастом и временем года. Более того, животные перемещаются. 
Характерная площадь местообитания составляет для моллюсков всего 1—100 м?, для мышевидных 
грызунов - 0,01-—1 га, для копытных — несколько квадратных километров. Отдельные особи, одна- 
ко, мигрируют соответственно на сотни метров, на километры и на десятки километров. Птицы 
мигрируют на сотни и тысячи километров. 

Пробоотбор животных в Зоне отчуждения не делали по регулярной сетке. Сборы моллюсков 
перекрывают только главные водоемы и водотоки. Сборы рыб, хотя и многочисленные, сделаны на 
немногих стационарах и в рыбохозяйственных прудах. Сборы наземных позвоночных оставляют мно- 
го «белых пятен» на планшете карты. Поэтому непрерывные карты загрязнения можно строить только 
путем экстраполяции закономерностей перехода радионуклидов из внешней среды в организм 
животных, если изменчивость и пространственное постоянство животных в принципе позволяют 
предсказать загрязнение образца из произвольной точки по известным образцам из соседних точек. 

Поставим задачи обработки радиометрических данных о животных для последующего картогра- 
фического отображения: 

1) уменьшить дисперсию данных за счет установления закономерных трендов накопления 
радионуклидов: возрастного, сезонного, многолетнего — и свести массив измерений для данного 
вида к стандартной дате, стандартному возрасту; 

2) уменьшить дисперсию данных за счет установления постоянных отношений между накопле- 
нием радионуклидов животными различных видов и свести массив измерений для группы видов к 
одному, принятому за стандарт; 

3) установить коэффициенты перехода и накопления, получив возможность предсказывать 
загрязнение животных по картам загрязнения среды и, обратно, судить о загрязнении почвы или 
воды по загрязнению индикаторных видов; 

4) установить статистическую зависимость между загрязнением животных, добытых в соседних 
или отдаленных точках; 

5) оценить надежность предсказания ожидаемого уровня загрязнения на основании 
статистической характеристики случайной изменчивости загрязнения у животных данного вида. 


Источники данных 


Границы планшета карты 85 х 65 км соответствуют стандарту (Атлас..., 1996). На карту нано- 
сятся данные о бета-активности раковин моллюсков как о показателе загрязнения г (Егапіѕеуісћ еї 
а|., 1996) и о загрязнении рыб, птиц, мышевидных грызунов и копытных !Сѕ. Использованы 
опубликованные результаты исследований в Зоне отчуждения и безусловного (обязательного) отсе- 
ления (ниже — Зона отчуждения) и базы данных Института зоологии НАН Украины (В. А. Гайчен- 
ко, Л. И. Францевич), Института зоологии АН Беларуси (А. И. Воронович, В. В. Ермолаев, П. Г. Коз- 
ло, С. В. Кучмель), Института гидробиологии НАН Украины (М. И. Кузьменко, Е. Н. Волкова, 
И. В. Паньков), Центра экологических проблем атомной энергетики Украины НЦ «Институт 
ядерных исследований» НАН Украины (О. Л. Зарубин), Института рыбного хозяйства УААН 
(И. А. Курочкин), лабораторий радиоэкологии и радиобиоэкологии леса Старопетровской лесной 
опытной станции по селекции и повышению продуктивности лесов, Полесской агролесомелиоратив- 
ной научно-опытной станции УкрНИИЛХ (В. П. Краснов, 3. М. Шелест). Обработка первичных 
данных выполнена автором статьи. 

Сопоставимые данные по загрязнению почвы в местообитаниях копытных были получены в 
ряде случаев в результате гамма-спектрометрии образцов почвы, отобранных в месте отстрела. Кроме 
того, использованы данные (Лихтарев и др., 1995; Радиационная..., 1995; Гриценко та 1н., 1998; 
Нагорський та ін., 1996; Кашпаров, 1998), база данных по Киевской обл. Украинского центра радиа- 
ционных и гидрометеорологических наблюдений. Данные о концентрации °5г в водоемах Зоны 
отчуждения заимствованы из: Войцехович та ін., 1996; Войцехович и др., 1997; Панасевич та ін., 
1998; Афонін та ін., 1998. 

Хотя факторы изменчивости рассматриваются порознь, но для их анализа мы часто привлекаем 
не первичные данные, а данные, стандартизованные после устранения тренда. Для стандартизации 
мы привлекаем следующие зависимости: 1) коэффициент накопления А„, или безразмерное отно- 
шение между концентрацией радионуклида в биологическом материале и воде; 2) коэффициент 
перехода / (ігап5їег Гасіог), или отношение между удельной активностью радионуклида в 
биологическом материале (кБк/кг) и поверхностной активностью радионуклида в 10-сантиметровом 
слое почвы (кБк/м?). Единицы ў (м?/кг) наглядно показывают, какая площадка почвы содержит в 
среднем столько же радионуклида, сколько в одном килограмме биологического материала. 

Укажем латинские названия исследованных видов животных: 


Стандартизация данных... 101 


Моллюски: большой прудовик — Гутпаеа ѕіарпаііѕ І іппаеоѕ. 

Рыбы: щука — Ебох Іисіцз Ілппаеџцѕ; густера — Вйсса Бјоеғкпа І лппаєця; чехонь — Ре/есиѕ сиЇїғаіиѕ 
|Лппаецз$; карп — Сургіпи5 сагріо Ілппаеџѕ; сом европейский — .5//иғиѕ ат Ілппаєця; карликовый 
сомик — Атіиғиѕ пебиіозиз е Ѕиеиг; судак — Гисюретса шсторегса ГАппаеиз$; толстолобик белый — 
Нурорћіћа/ітісћіпуѕ тоших Маіепсіеппеѕ; толстолобик пестрый — Алізпіситуз поз Кісһага. 

Птицы: тетерев — Гугигиу іеіғіх іппаеџиѕ; лысуха — Рийса ага І ппаец5; чайка озерная — Гагиз 
підїбипацз Ілппаеџѕ; кряква — Апаѕ ріаїугпупспоз Ііппаеиѕ; чирок-трескунок - А. диетдиеаша ГАппаецз, 
большая синица — Рағиѕ та)ок іппаеиѕ. 

Грызуны: полевка обыкновенная — Мстому агуа/іѕ РаПаѕ; полевка-зкономка — М. оесопотиѕ 
РаПаѕ; полевка рыжая лесная — С/ейптопотуз сіагеоіиѕ Ѕсһгебег; мышь полевая — Ароаетиз азтатих РаПаз. 

Копытные: кабан дикий — ,$и5 5сго/а І іппаеџѕ; лось европейский — А/сеѕ а/сеѕ іппаецѕ; олень 
благородный — Сегуиз е/арһиѕ Ііппаеџцѕ; косуля европейская — Саргео[иу саргеои5 іппаеиѕ. 


Возрастная изменчивость 


Известно, что зависимость загрязнения '!37С$ от возраста (или от массы те- 
ла) прямая у щуки (Шевцова, Воронович, 1995) и чехони (Насвит и др., 1997), 
но отсутствует у густеры (Зарубин, 1997). Существенные зависимости могут быть 
использованы для сведения данных к одному возрастному стандарту. Возрастные 
зависимости для копытных, грызунов и моллюсков не установлены. 


Сезонная изменчивость 


Накопление !Сѕ у косули и кабана связано с составом кормов в разные 
сезоны. Для обработки использованы не первичные значения радиоактивности 
мышечной ткани, а коэффициенты перехода {относительно загрязнения почвы, 
после нормализации путем преобразования 4 = ү. Годовая цикличность 4 ап- 
проксимируется первой гармоникой тригонометрического ряда: 


4 (р) = О + азшр + р соѕ р, (1) 

где О - среднее значение 4; р в градусах - фаза даты добычи (последняя 
выражена в днях Д от начала года), р = 270/365. Значения коэффициентов подо- 
браны методом наименьших квадратов и составляют: 

для кабана О = 1,523, а = 0,284, Б = 0,422; 

для косули О = 1,856, а = -0,764, В = 0,354. 

Первая гармоника исчерпывает соответственно у кабана 9,8%, у косули 
23,1% общей дисперсии 4. Уравнения сезонной изменчивости накопления отно- 
сительно накопления А, на стандартную дату 1 октября, где время / выражено в 
днях с начала года: 


для кабана А (1) = А, (1,235 + 0,230 ѕіп + 0,342 соѕ г); (2) 
для косули А (0) = А (0,710 - 0,292 ѕіп Е + 0,135 соѕ 7). (3) 


Результаты аппроксимации показаны на рисунке 1. Максимум (минимум) 
накопления приходится у кабана на 05.02 (05.08), а у косули на 20.10 (25.04). 


Рис. 1. Сезонные изменения коэффициента перехода 
137С$ у косули. График в полярных координатах пред- 
ставляет годовой цикл, месяцы указаны по секторам. 
Линии уровня 1—3—10—30 дм?/кг (преобразование ү). 
Кардиоида — лучшее приближение цикла, полученное 
методом наименьших квадратов для уравнения 
(уравнение 1). 


Ею. 1. Зеазопа| їгепа ої е їгапѕѓег стог ої С іп гое 
деег ( Сеғуиѕ еіарпиз). Ріої іп роіаг соогаіпаќѓеѕ гергеѕепіѕ 
Фе аппиаі сусіе, ѕесіогѕ питбегеа Бу топіћѕ. Геуе! Ппез 
Юг аге 0.01—0.03—0.1— 0.3 т2/ке ((гапѕѓогта(іоп (1/)05 ). 
Во! сагаіоіа ѕһоуѕ (һе Беѕі бі г е аппиа! сусіє, 
дегіуеа гот Ме їгівопотеїгіса! гом (едиаііоп гом 1). 
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Рис. 2. Распределение удельной бета-активности (кБк/кг) в раковинах большого прудовика (с. Копачи, 
база данных ИЗ НАНУ). п - количество образцов, т — среднее, У - коэффициент вариации (т/). 


Е. 2. Різітібийоп ої ће ѕресійс Беа-асНуйу (КВа/К») іп (ће ѕһе11 ої "Бе ропа зпаЙ Гутиаеа 5азпай5 
(затре Нот а стеек іп Корасфі, 5 кт гот Ше СҺМ№РР, даѓабаѕе ої 17, МАЗИ). п 15 питбег ої затріе8, 
т 15 Ше теап, И і5 ће соеїйсіепі ої уагіайоп (т/5). 


Рис. 3. Распределение скорости счета бета-частиц (имп./с) от черепов полевки обыкновенной 
(с. Копачи, 1987 г., данные ИЗ НАНУ). 


Во. 3. Юіѕігібибоп ої ће соипі гаѓе ої Бека рагіісІеѕ (ршѕеѕ/ѕ) іп ће агу 5К0115 ої "Бе соттоп уоіе Місғоѓиѕ 
ағуаііѕ (Корасһі, 1987, аѓабаѕе ої 17. МАЗИ). 


Рис. 4. Распределение удельной активности !37С$ (кБк/кг) в мышечной ткани кабана (Днепровско- 
Тетеревское ГОХ, база данных СПЛОС). 


Но. 4. Юіѕігібийоп ої (ће зресійс асііуїку ої З7С8 (КВа/Кке) іп (ће тизе (іѕѕие ої (ће Боаг 5и5 зсго/а. Тһе 
Ошщерег-Таетеу һипііпе Бизбапагу (дӢаѓабаѕе ої ће Ѕїагореїгоуѓѕі Еогеѕігу Кеѕеагсһ Ѕѓайоп). 
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Аналогичный осенний пик был отмечен у косуль в Баден-Вюрттемберге в 1987— 
1989 гг. по материалам отстрелов около 2400 животных (\№оГ 1990). В отличие 
от наблюдений в Чернобыле, в месте аварии на атомной станции на р. Саванна 
(США, штат Южная Каролина) годовой цикл у кабана был обратный: загрязне- 
ние летом было втрое больше чем зимой (ЗиШие е{ а|., 1986). Несходство 
можно объяснить разницей в широте и географической зональности, в климате, 
рельефе местности и, следовательно, в сезонном характере питания животных в 
Чернобыле и субтропической Южной Каролине. 


Внутривидовая изменчивость 


Измеряется в однородных выборках из одного местообитания, собранных в 
одно время или приведенных к стандартной дате, стандартному возрасту после 
устранения временных трендов. Примеры распределения индивидуальных уров- 
ней радиоактивности показаны на рисунках 2—4. Распределения имеют правую 
асимметрию. Форма распределения с параметрами т (среднее), с (среднеквадра- 
тическое отклонение), 5, (ошибка среднего) зависит от коэффициента вариации 
У = 100% • с/т. Хорошим приближением к реальным распределениям является 
логарифмически нормальное. 

Коэффициент вариации бета-активности раковин у пресноводных моллюс- 
ков из Зоны отчуждения составляет 10—30%, у сухопутных моллюсков 20—60% 
(Францевич и др., 1995). Так как любая проба содержит несколько экземпляров 
одного вида, ошибка измерения уменьшается за счет усреднения. 

Коэффициент вариации содержания !7С$ для густеры из водоема-охлади- 
теля (1990) составил 10—13% (Зарубин, 1997), для белого и пестрого толстолобов 
оттуда же (1987—1988) — 59% (база данных ИЯИ). В пойменных озерах (Смер- 
жов, Персток, Плоское) коэффициент вариации для мирных рыб в среднем 21% 
(пределы 3—56%), для хищных 46% (22—78%) (база данных ИЗ АНБ). 

В пределах планшета карты добыто около 200 образцов птиц, но сборы не 
были систематичными. Выбрано всего пять однородных групп (табл. 1). Чем боль- 
ше объем выборки, тем больше коэффициент вариации. Очевидно, верхние зна- 
чения (120—180%) ближе к истинным для больших популяций. 

Вариабельность грызунов установлена по достаточно представительным про- 
бам из одного места сбора (22—75 экз.) в пределах 50—180%. 

Коэффициент вариации содержания !37С$ у копытных составил после устра- 
нения временных трендов для кабана 134%, для косули 105%. Для ограниченных 
выборок лосей (13 экз.) и благородных оленей (54 экз.) из Днепровско-Тетерев- 
ского охотхозяйства коэффициенты вариации составили соответственно 127 и 
137%. Коэффициент вариации загрязнения косуль в Баден-Вюрттемберге по 
сборам летом или осенью (раздельно) в пределах одного лесничества находился 
в пределах 60—140% (среднее 86%) (МоІҒ, 1990). 

Коэффициент вариации минимален у оседлых животных, обитающих в 
однородной среде: у пресноводных моллюсков, у озерных рыб. Коэффициент тем 


Таблица 1. Внутривидовая изменчивость удельной активности !3?С5 (А„, кБк/кг) в мышечной ткани птиц 
Таріе 1. Іпігаѕресійс уагіарійсу ої 137С$ зрес с асйуНу (А,, КВа/Ке) іп Ше шизе И55ие ої Біга5 


Вид Место сбора Год А, среднее У, 96 п 
Кряква, чирки | ближняя зона 1987 5,726 179 28 
Кряква, чирки левый берег Припяти 1989--1990 4,875 54 9 
Большая синица Рьжий лес 1992 324,3 85 
Озерная чайка охотхозяйство «Вершина» 1986—1987 0,491 120 17 


Тетерев Машев 1990 98,01 34 6 
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о Рис. 5. Номограмма характеристик логнормаль- 

ного распределения: А/т — активность, отнесен- 

99 ная к среднему; О — квантиль распределения, %; 
У — коэффициент вариации, %; а — шкалы для 
0 < И< 100; Б — шкалы для 100 < И < 200. Пунк- 

95 тир — пример нахождения квантиля для И = 80%, 

90 А = 2т. 

80 Еів. 5. Моторгат оЁ ѕіабіѕісѕ ої (ће Іовпогта] 

70 аіѕігірибоп: А/т — асіїміїу, ге іо ће теап; О - 

60 ће диапійе ої ће аіѕігібийоп, іп %; И - сое сіепі 

50 оѓ уагіайоп, %; а — зса1ез ог 0 < И < 100; 6 - ѕсаІеѕ 

40 Юг 100 < И < 200. Раѕћеа Ппе (- -) - ап ехатріе 

30 ої ром ќо йпа ће дџапіе Гог И = 80%, А = 2т. 

20 

10 
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больше, чем больше животное склонно к миграциям. Погрешность прогноза 
уровня загрязнения отдельной особи не лучше, чем указанные значения коэффи- 
циентов вариации. С помощью номограммы (рис. 5) любой практический работ- 
ник, охотник или рыбак может, не прибегая к компьютеру, установить соотно- 
шение образцов, превышающих данное значение загрязнения, при известных зна- 
чениях среднего и коэффициента вариации. 


Межвидовая изменчивость 


Постоянные соотношения накопления %5г моллюсками разных родов рас- 
считаны для значений удельной бета-активности раковин в парных сочетаниях 
моллюсков из разных родов, собранных в одном месте. После логарифмического 
преобразования данных подбирались поправки, минимизирующие квадрат разни- 
цы между видами, одновременно для всех десяти родов. Один вид - большой 
прудовик, ГЃутпаеа 5іаєпаїіз — был принят за стандарт сравнения. Антилогарифм 
поправки для данного рода - это отношение накопления к стандарту (табл. 2). 
Без потери точности эти значения можно округлить до первой значащей цифры. 

Для сборов в пруде-охладителе были установлены такие соотношения между 
удельной активностью (А, ) '37С$ у отдельных видов рыб (Казаков и др., 1994): 


Таблица 2. Поправочные коэффициенты приведения бета-активности раковин моллюсков к стандарт- 
ному виду — большому прудовику, Гутпаеа ѕїарпаііѕ (база данных ИЗ НАНУ, 325 пар сравнения) 
Таре 2. Соггесіїоп диобіепія Гог гедисйоп ої Ше ѕресійс тоПизк зНе| Беёа-асйуйу о ће 5капдага, ће ропа 
5паї! (йаќараѕе ої 17, МАЗИ), 325 сопітопќей раігѕ ої ѕресіеѕ 


Род Коэффициент приведения Относительное накопление 
Тутпаеа 1,0000 1,0000 
РІапоғђағіиѕ 1,0101 0,9900 
Иуірағиѕ 0,9950 1,0050 
Апоаопіа 0,6281 1,5921 
Ото 1,0070 0,9930 
Фуеіззепа 0,5593 1,7879 
Сераеа 0,1572 6,3613 
Втадуфаепа 0,1534 6,5189 
Уиссіпеа 0,2783 3,5932 


Нейх 0,0436 22,936 
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карп 1,00 сом канальный 1,07 
белый толстолобик 0,91 сом европейский 3,20 
пестрый толстолобик 1,07 судак 3,04 

щука 5,40 


В сборах из других водоемов коэффициенты несколько отличаются от цити- 
рованных значений, но везде хищные рыбы (в период хронического поступле- 
ния радионуклида) загрязнены больше, чем мирные. Сопоставляя усредненные 
измерения на мирных и хищных рыбах, собранных в одно время в одном место- 
обитании (рис. 6), находим: 


ЇВ А, щн. = 0,9581 18 А, „ьн + 0,5134 (п = 28, ғ = 0,9606), (4) 


где п — количество измерений, г - коэффициент линейной корреляции. Округ- 
ленно среднее отношение накопления !Сѕ хищными и мирными рыбами состав- 
ляет 3,3. 

На упомянутом выше полигоне на р. Саванна коэффициенты накопления 
были выше, а соотношения между экологическими группами иные, чем в Черно- 
быле: 8500 у хищных рыб, 11 000 у планктонофагов, 6400 у всеядных рыб (Мћіс- 
Юег её а!., 1990). 

На фоне разброса загрязнения птиц не удалось найти систематической раз- 
ницы в загрязнении даже для таких массовых видов, как кряква и лысуха. 

Отношение накопления !Сѕ у кабана и косули (соответственно 195 и 136 осо- 
бей) установили сравнением коэффициентов перехода после устранения сезонного 
тренда. Отношение накопления на 01.10 составляет 0,2375, обратное отношение 
4,2. Для вычисления отношения накопления кабан/лось (1,53) и кабан/олень (1,18) 
использовали данные по Днепровско-Тетеревскому охотхозяйству. 


Коэффициенты перехода и накопления 


Распределение коэффициента перехода резко асимметричное; оно нормали- 
зуется после логарифмического преобразования. Большая вариабельность значе- 
ний ў зависит не только от изменчивости биохимических свойств животных, но 
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Рис. 6. Сопоставление удельной активности '7С$ (кБк/кг) в мышечной ткани мирных и хищных рыб 
из одного водоема (Шевцова, Воронович, 1995; Насвит и др., 1997; Вовк, 1996; база данных ИЗ 
АНБ). Данные в логарифмической шкале. 


Не. 6. Сотрагіѕоп ої ѕресійс асііуіќу ої '7Сз (КВа/Ке) іп ће тизе 115зие об ПитЫе (отпіуогоиѕ) апа 
ргедаїогу Из соПесѓеа іп һе зате уаег бойу (Шевцова, Воронович, 1995; Насвит и др., 1997; Вовк, 
1996; даѓабаѕе ої 17, МиззК). аа аге ргеѕепѓеа іп Іовагііһтіс зсае. 
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Рис. 7. Зависимость бета-активности раковин моллюсков А, (Бк/кг) (ориг.) от загрязнения %8г 
водных объектов А, (Бк/л), Зоны отчуждения; г — коэффициент корреляции. 


Но. 7. Рігесі гайоп беѓмееп %5$г сопсепігайоп іп (ће маїег А, (В4/1) іп (ће уаег бойіеѕ ої (ће Ехсазюп 
Гопе апа зресійс Беа-асНуйу ої тоПизК ѕһеѕ А, (КВа/Ка); г 15 согтейаНоп диойепі. 


и от неопределенности соотнесения места добычи животного с местом, где про- 
веден отбор образцов почвы. Значения бета-активности раковин моллюсков 
сопоставлены с опубликованными данными по загрязнению °%5г водных объе- 
ктов Зоны отчуждения. Данные 65 проб, приведенных к стандартному виду и 
дате, усреднены по 25 водоемам (рис. 7). Уравнение зависимости А,, Бк /кг, от 
А, Бк/л: 

16 А, = 1,1696 Ів А, + 3,5772 (и = 25, ғ = 0,9685), (5) 
где свободный член — логарифм коэффициента накопления, равного 3777. От- 
сюда можно вычислить, что для норматива на воду 2 Бк/л по г (Державні..., 
1997) критическая бета-активность раковин стандартного рода — прудовика — 
составляет 7,5 кБк/кг. 

Связь загрязнения местных моллюсков с известными данными по загрязне- 
нию почвы 95г (Лихтарев и др., 1995; Кашпаров, 1998; база данных Украинско- 
го центра радиационных и гидрометеорологических наблюдений) гораздо менее 
определенная. При использовании базы данных (Кашпаров, 1998), собранных по 
единой методике, или всех доступных данных получены одинаковые зависимости 
соответственно: 


16 А, = 0,9900 1 А, —0,9207 (и = 69, "= 0,7827); (6) 
16 А, = 0,9601 1 А, — 0,9111 (п = 105, г = 0,7579). (7) 
Свободный член — логарифм коэффициента перехода, равный для 


уравнения (6) 12 дм?/кг. Ранее для сборов моллюсков из Киевской области был 
получен средний по области коэффициент перехода для бета-активности рако- 
вин местных моллюсков (стандарт — прудовик) 11 дм?/кг (Францевич и др., 1995). 
Этот коэффициент, однако, менялся от северных районов области (20— 
30 дм?/кг) к южным (3—12 дм?/кг) обратно концентрации кальция в воде местных 
водоемов. В пределах Зоны отчуждения почвенные и гидрологические условия 
меняются мало. 

Коэффициент перехода относительно загрязнения почвы 37С5 для сухопут- 
ных оседлых птиц (тетерев из Машева) составил 2,5 дм?/ кг - величина того же 
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Таблица 3. Коэффицент перехода Г (т + 5,) относительно загрязнения почвы '7С5 для сборов 
мышевидных грызунов из ближней зоны 

Таріе 3. Тгапѕѓег Гасіог Г (т + 5„) гайуе (0 Ше 50Й сощаштайоп у 137С$ Гог годепі затр!ез от пе 
с105е5{ 7опе агоипа һе СИХРР 


Вид ї дм?/кг п 
Полевка-зкономка — Місгоѓіиѕ оесопотиѕ 21,23 = 3,02 25 
Мышь полевая — Ародетиз аетаних 5,56 0,73 133 
Рыжая лесная полевка — С/ейгуопоту5 &агеом5 14,97 + 3,34 64 


порядка, что и для копытных. Для уток из ближней зоны / может быть оценен 
в 0,38 дм?/кг на 1987 г., однако значение этой величины в настоящее время неиз- 
вестно, так как пробоотбора птиц в последние годы в Зоне отчуждения никто не 
проводил. Коэффициенты перехода для '7Сз для мышевидных грызунов даны в 
таблице 3. Активность '7С$ измерена в целых тушках грызунов. 

Связь между загрязнением почвы !Сѕ и загрязнением мышечной ткани 
копытных показана на рисунке 8 и в таблице 4. Кроме логарифмического пре- 
образования данных, для нормализации и графического отображения оказалось 
удобным также преобразование д = ў. Значения // (дм?/кг) до стандартизации 
показаны в таблице 4. При раздельной обработке украинских и белорусских сбо- 
ров оказалось, что коэффициенты для кабанов, полученные на обширном мате- 
риале, одинаковы (1,48 и 1,58; # = 0,54), мало различаются для косуль (1,94 и 
1,28; 1 = 2,97), но значимо различаются для лосей (1,02 и 0,56; ѓ = 11,6). Этому 
различию пока нет рационального пояснения. 


10000 Ат Рис. 8. Связь между загрязнением А, (кБк/м?) 
6 почвы '7С$ и загрязнением А, (кБк/кг) мы- 
1000 ө шечной ткани копытных (после устранения 
Я 8 сезонного тренда и приведения данных к 
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Таблица 4. Значения функций (т + 5,) от коэффициента перехода для '37С$, у копытных 
ТаЬ1е 4. Уаїше5 (теап т + 5„) ої Гипсбіоп5 ої ёһе (гапяїег ѓасіог 1/ ѓог !37С5, аш?/К?, іп ипешаќеѕ 


Выд Функция от коэффициента перехода при преобразовании, дм?/кг М 
х= у % х= 15 ў ха 
Лось 1,0286 = 1,5526 151 —0,3332 + 0,5522 0,8365 + 0,5734 50 
Кабан 3,6763 + 5,3555 146 0,0944 + 0,7206 1,5232 + 1,1646 195 
Олень 2,2386 + 3,0795 138 0,1489 + 0,3925 1,3220 + 0,7007 54 
Косуля 4,8993 + 5,9139 121 0,3290 + 0,6383 1,8562 = 1,2057 136 


Стандарт 2,0482 + 3,8218 186 0,0482 + 0,6129 1,1835 + 0,8047 436 
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Рис. 9. Изменение загрязнения '7С$ в мирных (1, 3, 5) и хищных (2, 4, 6) рыбах р. Припять (1, 2), 
Киевского водохранилища (3, 4) и пруда-охладителя (5, 6) в 1986—1995 гг. (Насвит и др., 1997; Вовк 
та ін., 1996). 


Рр. 9. Тгепа оѓ З"С58 сопіатіпа(іоп дигіпе 1989—1995 іп (ће һҺитЫе (1, 3, 5) апа ргейаќогу (2, 4, 6) БП 
тот фе пуег Ргіруаї (1, 2), Куіу Кезегуо (3, 4) апа ће СҺМРР содіїпе ропа (5, 6). (Егот Насвит и 
др., 1997; Вовк та ін., 1996). 


После стандартизации данных по всем копытным к виду и сезону и сведе- 
ния всех измерений в единый список получены следующие прямые и обратные 
уравнения регрессии (п = 436): 


ів А, = 0,6820 № А, - 1,2064 ("= 0,6909); (8) 
16 А, = 0,6999 18 А, + 2,2279 ("= 0,6909); (9) 
їв = —0,3180 Ів А, + 0,7935 ("= 0,4070). (10) 


Значения функций коэффициента перехода для стандартного копытного 
(кабан, 01.10) и различным образом преобразованных данных показаны также в 
таблице 4. Значения для линейного, логарифмического и квадратного корня 
преобразований составляют 2,05; 0,90 и 1,40 дм?/кг. Сходное по порядку вели- 
чины значение // - 3,6 дм?/кг исследователи получили для копытных на загряз- 
ненных территориях Беларуси в 1992—1993 гг. (КетиезБег» ей а1., 1996). Значение 
у около 2 дм?/кг можно рассчитать по данным (УМ/0ЇЇ, 1990) для косуль в Баден- 
Вюрттемберге. 


Многолетний тренд 


По 32 парам сборов моллюсков в одном и том же месте в 1988—1994 гг. и 
1997 г. установлены постоянная времени спада бета-активности раковин 
16,8 + 4,3 года (т + 5,) и период полуочищения — 11,6 лет. По трем стациона- 
рам (пруд-охладитель, Припять, Киевское водохранилище) установлен много- 
летний тренд загрязнения рыб '37С$ (Насвит и др., 1997; Вовк та ін., 1996). По- 
стоянная времени спада установлена в 4,022 года, период полуспада 2,8 лет 
(рис. 9). По литературным данным (Вгіѕбіп ег а!., 2002), период полуспада загряз- 
нения рыбы /еротіѕ тасғоһіғиѕ З"С8 находился между 3 и 7 годами. 

Для копытных (после стандартизации удельной активности !7С$ к кабану на 
01.10) корреляция величины 4 = У с годом добычи составила —0,0824. Много- 
летний тренд загрязнения для этой группы животных отсутствует. 


Стандартизация данньх... 109 


Стандарты накопления 


Из большого количества видов животных, обитающих в Зоне отчуждения, 
мы отобрали несколько видов, представительных для отрядов или классов 
животных, которые входят в пишевой рацион человека как объекты охоты и 
рыболовства и/или являются удобными индикаторами загрязнения. Эти виды — 
кабан, серая полевка, кряква, карп, большой прудовик. Сборы птиц в Зоне от- 
чуждения охватывали первые послеаварийные годы, а сборы мышевидных гры- 
зунов проводили в немногих стационарах. 

Авторам цитированых статей и автору настоящей статьи удалось установить 
сезонные тренды радиоактивного загрязнения для кабана и косули, многолетний 
тренд для моллюсков и рыб, стабильные отношения накопления у разных видов, 
родов, экологических групп у моллюсков, рыб, копытных. На основании назван- 
ных коэффициентов и функциональных зависимостей данные измерений были 
пересчитаны на стандартный вид и дату І октября 1996 г. 

Для моллюсков приведение к стандартному виду уменьшает разброс изме- 
рений в многовидовых сборах: после введения поправок начальная дисперсия 
логарифмически преобразованных данных уменьшается на 51%. 

Коэффициент вариации // после стандартизации уменьшился для линейных 
данных у кабана со 146% до 134%, у косули со 121% до 105%. Среднеквадрати- 
ческое отклонение 1 А, первичных данных составляет 0,8475 лог. ед., после 
стандартизации к виду и сезону оно составляет 0,7743, а среднеквадратическое 
отклонение 12 // до стандартизации составляет 0,6753, после стандартизации — 
0,6129. Таким образом, операция стандартизации устраняет 17% общей диспер- 
сии первичных данных, операция вычисления коэффициента перехода — 37%, 
обе вместе они уменьшают дисперсию данных на 48%. 


Пространственные закономерности 
радиоактивного загрязнения 


Зависимость А,, кБк/кг, от расстояния Џ, км, от ЧАЭС для моллюсков из 
местных водоемов и для копытных описывается соответственно уравнениями: 

1 А, = —1,8815 Ів р + 3,2787 (л = -0,6102, и = 165); (11) 

16 А, = -0,8322 Ів р + 1,8849 (г = -0,2927, п = 410). (12) 

Низкий коэффициент пропорциональности для копытных объясняется 
большой изменчивостью данных. У менее вариабельных моллюсков коэффици- 
ент пропорциональности приближается к —2 (загрязнение обратно пропорцио- 
нально квадрату расстояния, рис. 10). 

В действительности радиоактивное загрязнение местности еще и пятнистое. 
Необходимо установить, можно ли по дискретным пробам животных строить 
содержательную непрерывную карту загрязнения, т. е. прогнозировать загрязне- 
ние объекта в произвольной точке карты по данным пробоотбора в соседних или 
более отдаленных точках. 

Надо рассчитать сходство измерений в соседних точках в зависимости от 
расстояния между ними (произведение значений измерений или корреляцион- 
ную функцию). Предварительно данные о загрязнении после логарифмического 
преобразования нормализуются: вычитается среднее значение, остаток делится 
на среднеквадратическое отклонение. Корреляционная функция — это сглажен- 
ное значение произведения пары измерений в зависимости от расстояния между 
точками пробоотбора. 

Пример пространственного анализа данных о загрязнении копытных в Дне- 
провско-Тетеревском охотхозяйстве (рис. 11) показывает, что корреляционная 
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Рис. 10. Обратная зависимость между расстоянием от ЧАЭС (р, км) и удельной бета-активностью 
раковин моллюсков (А,, кБк/кг), стандартизованной по прудовику. 


Ее. 10. [шуетзе геіаіопѕһір Бебмеєп {ће аіѕіапсе Нот Ше СЬМРР (О, кт) апа ѕресійс Беа-асйуйу ої 
тойизк ѕһе115 (А,, КВа/Ке.), ѕіапаагаіхеа го пе ропа зпай. 


Рис. 11. Сходство между значениями загрязнения '7С$ проб копытных в зависимости от расстояния 
ДО между точками пробоотбора в Днепровско-Тетеревском охотхозяйстве, оцененное как 
произведение Р нормализованных значений загрязнения в двух точках. Произведения парных 
значений 122 пробы, 7323 пары сравнения. 


Ее. 11. 5ітіагіїу Бебуееп пеірһбог ипешаѓе затр!ез ут Фе Ощерег-Теегеу рипіїпе Һиѕбапагу, 
ехргеззе4 аз ће ргодисі А ої погта|2е уашез ої сопгатіпайоп мії '7С$ іп а раїг ої ѕатріеѕ уегѕиѕ 
дізгапсе Бебмееп Бо ѕатрііпе ѕіќеѕ (122 затр|ез, 7323 сопйоще4 раігѕ). 


Рис. 12. Корреляционная функция К бета-активности раковин моллюсков в зависимости от 
расстояния между точками пробоотбора в окрестностях Зоны отчуждения (165 проб, 13 435 пар 
сравнения). 


Ее. 12. Соггеіайїоп Гипсбіоп В Бебмееп {ће ѕресійс зВеЙ Беїа-асіуіу іп раігѕ ої тоШазК ѕатріеѕ уегѕиѕ 
аіѕќапсе Бебмееп затрИи» ѕііеѕ іп усшез ої фе СҺМРР ($апдаг412е4 Чака, 165 затр!ез, 13 435 
сопітопіед раігѕ). 


Стандартизация данньх... 111 


функция может бьть близкой к нулю, независимо от расстояния. Для зтой 
территории нельзя выделить участки (пятна, зоны) с различным загрязнением и, 
следовательно, нет достаточных оснований строить непрерывную карту. Един- 
ственным разрешенным прогнозом загрязнения зверей является указание сред- 
него значения и его ошибки для всей территории. 

Оперируя с большей территорией и большим диапазоном загрязнения в пре- 
делах планшета карты, получаем статистически существенные значения корреляции 
для местных моллюсков (рис. 12) и копытных. Для сборов из одного и того же 
местообитания коэффициенты корреляции соответственно равны 0,89 и 0,59 (13 435 
и 83 845 пар сравнения). Если устранить тренд по расстоянию (уравнение 12) из 
измерений копытных, то значения корреляционной функции уменьшаются (для 
одного местообитания — до 0,44), однако ход функции не изменяется. Область 
положительной корреляции окружена периферией с отрицательной корреляцией 
относительно центра, т. е. «чистые» и «грязные» пятна чередуются. И для моллюс- 
ков, и для копытных радиус положительной автокорреляции около 30 км. 


Выводы 


Установлены сезонные тренды радиоактивного загрязнения для кабана и ко- 
сули, многолетний тренд для моллюсков и рыб, стабильные отношения накопле- 
ния у разных видов, родов, экологических групп у моллюсков, рыб, копытных. 

Максимум (минимум) накопления !7Сз приходится у кабана на 05.02 
(05.08), а у косули на 20.10 (25.04). 

Период полуочищения раковин моллюсков от загрязнения %8г составил 
11,6 лет, мышечной ткани рыб от загрязнения !37С$ 2,8 лет. 

На основании известных коэффициентов и функциональных зависимостей 
данные измерений были пересчитаны на стандартный вид (для моллюсков — 
большой прудовик, для рыб — карп, для копытных — кабан) и стандартную дату 
01.10.1996. Стандартизация данных вдвое уменьшает первоначальную дисперсию 
измерений в данном местообитании. 

Индивидуальный (внутривидовой ) коэффициент вариации (%) наименьший 
у оседлых гидробионтов — моллюсков (10—30) и рыб (10—60), больше у оседлых 
наземных моллюсков (20—60) и мышевидных грызунов (39—180) и максимален 
у способных к более или менее отдаленным миграциям копытных (120—150) и 
тем более птиц (120—180). Предсказание значений радиоактивного загрязнения 
животных с меньшей ошибкой нереально. 

Радиоактивное загрязнение животных убывает с увеличением расстояния от 
ЧАЭС, по-видимому, обратно пропорционально квадрату расстояния. Эта зави- 
симость зашумлена индивидуальной изменчивостью и, вероятно, пятнистым 
характером загрязнения местности. 

Радиоактивное загрязнение животных прямо пропорционально загрязнению 
местности или воды. Коэффициент накопления бета-активности (°5г/У) в ракови- 
нах моллюсков относительно концентрации %$г в воде составляет 3770, коэффи- 
циенты перехода '37С$ относительно загрязнения почвы составляют у наземных птиц 
2,5, у мышевидных грызунов 5,5—21, у стандартного кабана 2,05 дм?/кг. 

Обнаружена положительная корреляция загрязнения проб из соседних точек 
местности, тем большая, чем меньше расстояние между пробами. Для моллюс- 
ков и копытных радиус положительной автокорреляции 30 км. Этот радиус ха- 
рактеризует средний размер пятна загрязнения. 

На карту были вынесены две категории переменных: 

1. Карта фактов содержания !Сѕ в мышечной ткани стандартной рыбы, 
птиц, грызунов, стандартного копытного и бета-активности раковины стандарт- 
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ного моллюска как индикатор загрязнения %$г. Уровни загрязнения градирова- 
ны по логарифмической шкале (две градации на десятичный порядок). Второй 
уровень отражает допустимый уровень содержания С в рыбе и мясе копытных 
и бета-активность раковин, соответствующую допустимому уровню 9085г в воде 
(Державне..., 1997). Для Зоны отчуждения ЧАЭС эти карты были опубликованы 
в статьях (Францевич, Іщук, 2000; Францевич та ін., 1998); 

2. Непрерывная прогностическая карта загрязнения С стандартного каба- 
на в сухопутных ценозах и раковин стандартного прудовика %$г в водоемах. 
Принимая во внимание размер участка обитания у копытных, мы решили вы- 
полнить карту не в изолиниях, а ограничиться разбиением ее на квадраты 
4 х4км, в соответствии с координатной сеткой карты, принятой в Атласе (Ат- 
лас..., 1996). Для Зоны отчуждения карта опубликована в Бюлетені екологічного 
стану Зони відчуження (Францевич, Пцук, 2000); полная карта масштаба 
1: 200 000 в виде распечатки или электронной копии может быть заказана в 
ИнфоЦентре «Чернобыль» (Е-таїї: ѕсіепі@іс-сһегпобуі.Кіеу.џа.). Образец карты 
экспонируется в музее «Чернобыль» в Киеве. 
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